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Abstract 

Optical properties of a plasmonic nano-antenna made of two metallic nanospheroids (prolate or 
oblate) are investigated analytically in quasistatic approximation. It is shown that in clusters of 
two nanospheroids, three types of plasmonic modes can be present. Two of them can be effectively 
excited by a plane electromagnetic wave, while the third one can be effectively excited only by a 
nanolocalized light somxe (an atom, a molecule, a quantum dot) placed in the gap between the 
nanoparticles. Analytical expressions for absorption and scattering cross-sections, enhancement 
of a local field, and radiative decay rate of a dipole source placed near such a nano-antenna are 
presented and analyzed. 

PACS numbers: 73.20.Mf, 78.67.-n 
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I. INTRODUCTION 



Nowadays, quite a number of works are devoted to the study of optical properties of single 
nanoparticles and their clusters. Special attention is paid to metal nanoparticles with the 
help of which it is possible to enhance electric fields at the frequencies of localized plasmon 
resonances^"-. On the basis of this effect, a variety of possible applications are considered. 
The most developed is use of large local fields near a rough surface for increasing of surface- 
enhanced Raman scattering (SERS)^. Modification of fluorescence by means of nanoparticles 
of different shapes is a basis for creation of nanobiosensors^""-, nano-antennas^"— , devices 
for decoding of DNA structure^^, and etc. 

At the present time, optical properties of single metallic nanospheres and their clustersi^"— , 
single nanospheroids^"— , nanoellipsoids22,"i2^^ and some other nanobodies^>^>^ are studied 
well enough from analytical point of view. Many other nanoparticles' shapes and also 
nanoparticles clusters are investigated only numerically. Unfortunately, numerical simula- 
tions often do not allow understanding of physical nature of interesting and complicated 
phenomena in this area. That is why any analytical solutions will have principal importance. 

In this work, we present the results of an analytical study of optical properties of clusters 
of two metallic prolate or oblate spheroidal nanoparticles. Such clusters are investigated both 
experimentally and numerically and form a basis for various possible applications including 
nanosensors, nano-antennas, and plasmon waveguides^^"— . As far as we know, there is no 
analytical investigation of optical properties of two-nanospheroid clusters. That is caused by 
extremely cumbersome mathematical calculations connected with use of general spheroidal 
wave functions and general rotational-translational addition theorems^"—. However, in the 
case of nanoparticles which sizes are less than the wavelength, i.e. in a quasistatic case, 
one can neglect retardation effects. It allows us to find an analytical solution for a cluster 
of two spheroidal nanoparticles, placed in an arbitrary external field. The geometry of the 
considered problem is shown in Fig. [H For simplicity, we will consider that the cluster 
consists of two equal spheroids made of a material with dielectric constant e and placed in 
vacuum. Generalization of our approach to unequal spheroids is straightforward and can be 
easily done with the help of the general theorem presented in the Appendix. 

Main attention will be paid to a case of nearly touching and strongly interacting spheroids 
since it is the case that seems to be most interesting for applications because a substantial 
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enhancement of electric fields occurs there. The opposite case of weakly interacting spheroids 
can be easily treated with approximation of the spheroids by point dipoles with correspond- 
ing polarizabilities^. 

For illustration of the analytical results obtained, we will consider a case of two identical 
(prolate or oblate) nanospheroids made of silver—. We suppose that the largest size of the 
nanospheroid is equal to 30 nm, and the aspect ratio of the spheroid is taken to be equal to 



The rest of the paper is organized as follows. In the Section [TTl free plasmon oscillations 
of a two-nanospheroid cluster are investigated. The results of this section reveal underlying 
physics and are necessary for interpretation of results of the other sections. In the Section [ITU 
we will consider optical properties of the two-nanospheroid cluster placed in the field of a 
plane electromagnetic wave. Here, we will find absorption and scattering cross-sections and 
the factor of local field enhancement also. In the Section IIVI an object of examination is 
optical properties of the two-nanospheroid cluster placed in the field of a radiating atom 
or a molecule, and their decay rates will be calculated. In the Appendix, derivation of the 
translational addition theorem for spheroidal functions in a quasistatic case is presented. 
It is the theorem that allows one to find an analytical presentation of the solution for the 
two-nanospheroid cluster. 

II. PLASMON OSCILLATIONS IN A CLUSTER OF TWO NANOSPHEROIDS 

It is well-known that all optical properties of nanoparticles can be derived from their 
plasmonic spectra, that is, from related plasmon eigenvalues and eigenfunctions e,^, h,^, 
that are solution of the sourceless Maxwell equations^: 



where e = inside the nanoparticle and e = 1 outside it, w is frequency of electromagnetic 
oscillations and is the speed of light in vacuum. As a result, electric field in the presence 
of any nanoparticle can be presented in the following form^: 



0.6. 




(1) 
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where e (a;) describes the dependence of dielectric permittivity of nanoparticle's specific 
material on frequency u, Eq is the excitation field, and u is a vector index that defines 
specific plasmonic mode. From Eq. ([2]), it is possible to find any optical properties of 
a nanoparticle or a cluster of nanoparticles. So, to understand very complicated optical 
properties of a two-nanospheroid cluster, we should investigate plasmonic spectrum of this 
system at first. 

For studying plasmon oscillations and other optical properties of clusters of two nanospheroids, 
it is enough to solve the quasistatic equations: 

div (ee^) = 0, rotej^ = 0, (3) 
that can be reduced to solution of Laplace equations by substituting e^, = —Vipu- 



\s 



dn 



S 



Af/)^" = 0, inside the nanoparticle, 
Ay9°"* = 0, outside the nanoparticle, 
at the surface 
of the nanoparticle. 



in"'- 
Vu l5' 



(4) 



In Eq. (jl]), ip''^, (fl^^ are the potentials of plasmonic eigenfunctions inside and outside the 
nanoparticle correspondingly, and denotes normal derivative at the nanoparticles' 

surface S. The last equation in (j4]) provides continuity of normal components of electrical 
induction. Note that in this case there is no need to find magnetic fields for description of 
plasmonic oscillations. 

The systems of equations obtained in such a way, have nontrivial solutions only for some 
negative values of permittivity e^, defining frequency of plasmon oscillations^i^. In the case 
of Drude theory, e (u) = 1 — Upju"^, frequency of plasmon oscillations can be found from 
the expression 



pi} 



(5) 
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where oopi is bulk plasmon frequency of a metal from which the nanoparticles are made. 
Our approach allows us to investigate arbitrary spheroids, but for simplicity in the present 
section we examine equations for plasmon oscillations in a cluster of two identical metal 
nanospheroids. 

In our case of a two-nanospheroid cluster, we will look for a solution as follows. The total 
potential outside the spheroids will be the sum of their partial potentials (we will omit the 
mode index u further )^'^: 

,„OUt ,„OUt I ,„OUt /'f}\ 

while the potentials inside each nanospheroid shall be denoted as (j = 1, 2). To find 
^out ^in^ ^in^ naturally to use spheroidal coordinates. In the case of a prolate 
nanospheroid, the relation between Cartesian and spheroidal coordinates (1 < ^ < cxd, 
-1 < r/ < 1, < < 27r) is^: 



a: = /V(e-l)(l-^')cos0, 
y = /V(e-l)(l-r/2)sin0, 

z = Kv, (7) 

where / = \/ (? — (P' is a half of the focal distance in a prolate spheroid (a < c) which surface 
is set by the equation (x^ + y^) /a^ + 2;^/c^ = 1. 

In the case of an oblate spheroid (a > c), the relation between Cartesian and spheroidal 
coordinates (0 < < oo, — l^^^^l, O<0< 27r) has the following form^: 



X = /v/(e + l)(l-^')cos0, 
y = /V(e' + l)(l-^')sin0, 

^ = Kv, (8) 

where / = \/ — (? is a half of the focal distance in the oblate spheroid. Let us note that 
this expression can be obtained from Eq. (JTj) by substitution ^ and / — t- —if . Further, 
we will use this formal replacement since it is fundamental and allows us to find a solution 
for oblate spheroids if the solution for prolate spheroid geometry is known^"— . 
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A. Plasmon oscillations in a cluster of two identical prolate nanospheroids 



To find plasmonic spectra of a two-nanosplieroid cluster, it is naturally to use two local 
systems of spheroidal coordinates {^j, rjj, (pj, j = 1, 2) the origins oj of which are placed 
in the centers of corresponding nanospheroids (see Fig. [T]^a)). Coordinates (and all other 
values) related to the first or second nanospheroid will be denoted by the index "1" or "2" 
respectively. The potential inside the j-th nanospheroid can be presented in the form^ 
(j = 1, 2): 

oo n 

= Y.Y.Pn {Q Pn iv,) COS (m0,) + Sin (m0,)) , (9) 

n=0 m=0 

where [rj) is an associated Legendre function^ defined in the region — 1 < r/ < 1, {^) 
is an associated Legendre function^ defined in a complex plane with the branch cut from 
— oo to +1. The partial potential outside the j-th nanospheroid can be presented^ as {j = 1, 
2) 

oo n 

= J2T.Qn (0) Pn iVj) {C^JI COS (m0,) + Dlii sin (m0,)) , (10) 

n=0 m=0 

where {^) is an associated Legendre function of the second kind^ defined in a complex 
plane with the branch cut from — oo to +1. 

By construction, the potentials iQ and ffTOl) are solutions of the Laplace equation^^. So 
to find a solution of Eq. (jlj), one should use only the boundary conditions: 



in I , ^out I 



Q^OUt 

a=«o ~ ^ 

Q^OUt 



€i=«o 



6=Co 



(11) 



where = c/ \/ — a? = c/ f are local radial coordinates defining surfaces of the nanospheroids, 
e is permittivity of materials from which the nanoparticles are made. To reduce the bound- 
ary conditions ( iTTl) to a system of linear equations, we apply the translational addition 
theorem for wave functions of the prolate nanospheroid (see the Appendix). In the case of 
two identical coaxial nanospheroids, this theorem gives (j, s = 1, 2, j 7^ s, 0i = ^2) 



Qn (0) Pn iVj) = E ^S-n'" (^s) P^ iVs), (12) 
q=m 

where the functions Smlmn = S.mqmn if, f, I, 0) and Smqmn = Smqmn (/, /, 1-, 71") are defined in 
the Appendix. Applying the boundary conditions f|TT]) and the theorem ( !T2|) . one can obtain 
the following system of equations (n = 0, 1, 2, . . . ; m = 0, 1, 2, . . . , n): 



dP^ (^o) X TD-m (c \ dQn (^o) \ ^(1) 

+ - 1) PT (^o) ^^^^E 5i°L. (-1)"^^^ c^S 



dPn{ia)^m/^\ T-,m f ^ \ dQn iio) 



e- 



Qn (eo) - Pr (6) (-1)™+" 



1) Pn (eo) ^^^E ^i°L.C^S = 0- (13) 



q=m 



^(2) 

-'mnmq 



When deriving Eq. (fT3|) . we make use of the fact that for identical nanospheroids S„ 
(-1)'^'^'^ Smlmq (scc the Appendix) and take Smlmq = (-1)'"'*''^ Smnmg- The system of equa- 
tions for Dmn is identical to ( fT3l) and gives no additional information for plasmonic spectra 



of coaxial spheroids. So, we will not consider it further. 

As it results from symmetry of the considered cluster and the system (fT3!) . there are 
two independent types of solutions (plasmonic modes) with opposite parity. To select the 
first type of the modes (symmetric in ^ — )■ —z transformation), one should choose Cml. = 
(— 1)™ " Cmn in Eq. ( |T3|) . As a result, we shall obtain the following system of equations for 
the symmetric modes: 



dP^ (^o) ^m fc \ um (c \ (^o) \ ^(1) 

— ^^0 — ~ — d^o — / " 

+ - 1) Pn (eo) ^^^E Sl^LqCi^l = 0. (14) 



q=m 



To obtain the second type of the modes (antisymmetric in z ^ —z transformation), one 
should put Cmn = — (—1) Cmn in Eq. (fT3|) . As a result, we shall obtain the following 
system of equations for the antisymmetric modes: 
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( 



Qn (eo) - Pn (eo) 



dQ": (eo) 



) 




dp::- (eo) 

d^o 



oo 



ie - 1) PT (^o) 



Ec(o) Ml) ^ Q 



(15) 



q=m 



It is important to notice that separation of spectra into symmetric and antisymmetric 
plasmon modes is possible only in the case when there is a plane of symmetry. When m is 
even, antisymmetric modes have nonzero dipole moment and they are "bright" modes. On 
the contrary, symmetric modes have zero dipole moment and are "dark" modes when m is 
even. In the case of odd m, the "bright" and "dark" modes correspond to the symmetric 
and antisymmetric modes respectively. One can expect that the antisymmetric mode m = 
will have the largest polarizability and thus will be the "brightest" one for excitation of our 
cluster with a longitudinally (along the z axis) polarized plane wave. 

To study plasmon oscillations in clusters of two prolate spheroidal nanoparticles, we have 
solved the eigenvalue problems ( fT4l) and (|T5l) numerically. In Fig. [2], normalized plasmon 
frequency ou/upi of a cluster of two prolate nanospheroids (see Fig. [T](a) for the geome- 
try), corresponding to the first three plasmon modes, is shown as a function of normalized 
distances I /2c between the nanoparticles' centers. Eigenvalues e have been obtained as a 
nontrivial solution of the equations systems (HM and ( |T5|) in the case of an axis-symmetric 
problem (m = 0). Then, the found solutions have been substituted into Eq. (jS]) to obtain 
plasmon oscillations frequency. 

In Fig. 121 one can observe that plasmon frequencies of a cluster of two prolate nanospheroids 
tend to plasmon frequencies of a single nanospheroid (see Fig. W^c)) if the distances between 
the nanospheroids are large enough. When width of the gap between the nanospheroids 
tends to zero, solutions of the equations f|T^ and f|T5|) behave very differently. For sym- 
metric modes (Fig. EJ^a)), there are two branches: T-modes and M- modes. Modes of "T" 
type can be obtained by the method of hybridization of plasmon modes of a single prolate 
nanospheroid, analogously to the case of a cluster of two spherical nanoparticles^^. When 
width of the gap between the nanoparticles is decreasing to zero, normalized plasmonic 
frequencies of T-modes tend to various values in the range from to 1/ ^/2, on the analogy 
with a two-sphere cluster— i^. T-modes with higher indices (not shown for clarity) will 
concentrate near cu/cOpi = l/\/2. In Fig. [2t^a), one can also see that at very short distances 
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between the nanospheroids (//2c < 1.1), a new type of plasmonic modes (M-modes) ap- 
pears. M-modes are characterized by strong spatial locahzation in the gap between the 
nanoparticles. As a result, they can be effectively excited only by a strongly nonuniform 
electric field of the molecule or the quantum dot. Values of plasmonic frequencies of these 
modes lie in the range ujpi/\/2 < to < copi. As the gap width decreases to zero, plasmon 
frequency of M-modes tends to bulk plasmon frequency Upi. 

In Fig. El^b), nontrivial solutions of the equations system ( ITSll for the antisymmetric po- 
tential in an axial-symmetric case (m = 0), are shown. By analogy with a two-sphere cluster, 
we will call these modes L-modes (longitudinal) because they are "bright" only for longi- 
tudinal excitation. These modes can be described by the hybridization method of plasmon 
oscillations of single nanospheroids forming the considered cluster. As width of the gap 
between prolate nanospheroids decreases to zero, normalized plasmon frequencies of these 
modes tend to zero as it also takes place in the case of spherical nanoparticles^'^. Plasmonic 
frequencies of L-modes of higher orders (not shown) tend to Upi/\/2, and concentration of 
infinite number of L-modes occurs near this value. 

In Fig. [HI distribution of a surface charge of plasmonic modes of the lowest order in 
clusters of two identical prolate nanospheroids, is shown. It is seen in this figure that 
the T- and M-modes have symmetric distribution of the surface charge in contrast to the 
antisymmetric L=l mode. This behavior, of course, is in agreement with symmetry of the 
equations f|T4|) and f|T5|) . Another interesting feature is that the surface charge of T-modes is 
distributed over the surface of the all nanoparticles for any distances between them, while for 
M- and L-modes it is concentrated near the gap between the nanospheroids if the distance 
between them is sufficiently small. It is interesting to note also that the surface charge 
of M-modes is more concentrated in comparison with that of L-modes. Indeed, due to an 
electroneutrality requirement, the total surface charge on each nanospheroid should be equal 
to zero. Here, both positive and negative charges of M-modes are localized near the gap 
between the nanoparticles so that on the rest of nanoparticles the charge is almost equal 
to zero, as it is well seen in Fig. [3]^a). At the same time, in the case of L-modes for each 
of nanospheroids near the gap a charge of only one sign is concentrated, and a charge of 
the opposite sign is distributed with small magnitude over the remaining surface of the 
nanoparticles. Therefore, strictly speaking, the surface charge in an L-mode is distributed 
over the whole surface of the cluster nanoparticles though it is not clearly seen at small 
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distances between the nanoparticles (see Fig. [3]^a)). As the distance increases, the charge 
distribution changes in the cluster: it spreads over the nanoparticles' surface, tending in the 
limit to a distribution corresponding to single prolate nanospheroids (see Fig. [31(c)). 



B. Plasmon oscillations in a cluster of two identical oblate nanospheroids 

In this geometry, one should also use local systems of coordinates (^j, rjj, j = 1, 
2) which are connected with each nanospheroid and have origins Oj in their centers (see 
Fig. mb)). Now, the electric potential inside the j-th nanospheroid can be presented in the 
form (j = 1, 2): 

CXD n 

= E E ^n"^ (^^^•) iv,) {411 COS (m0,) + 5^) Sin (m0,)) , (16) 

71=0 m=0 

and the partial potential outside the j-th oblate nanospheroid will look like (j = 1, 2): 



oo n 

= E E {^Q Pn iVj) {C^dl COS (m0,) + D^^ sin (m0,)) . (17) 

The total potential outside the nanospheroids will be expressed by Eq. (E]). As boundary 
conditions for the potential, Eq. ( ITT]) , where = c/\/a^ — = c/ f, is used. In the case of 
oblate nanospheroids, the addition translation theorem (see the Appendix) has the following 
form (j, s = 1, 2, s)- 



Q- m p: m { -±±p; (>«.) p^ (..) 'f" , (18) 

sin(m0j) q=o p=o N^qinn sin \ 



where the functions M^]lnn = Mpqmn{-if,-ifJ,T''), N^lLn = Npqmn{-if,-'ifJ,TT) and 
Mpqlnn = Mpq^n {-if, -if, I, 0), A^'^gmn = Npqmn {-if, -if, I, 0) are defined in the Appendix. 
Now, substituting Eqs. ( fT6ll and (|T71) into Eq. ( fTTl) and making use of the addition translation 
theorem ( !T8|) . we shall obtain the following system of equations (n = 0, 1, 2, . . . ; m = 0, 1, 
2, n): 
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+ m ^^^f: E m!^1„ i-lf Cf^ = 0, 

+ ie-l) P- {^^o) ^^^f: E = 0. (19) 

By deriving Eq. ( IT^ . we take into account that Mmlpq = (— 1)""^^ Mmnpg (see the Ap- 
pendix) and denote Mmnpq = {— I)™' Mmlpq- The system of equations for Dmn is analogous 
to Eq. f|T9l) . and we will not analyze it here. 

Due to symmetry of a cluster of two identical oblate nanospheroids, there are two types 
of plasmon oscillation: symmetric and antisymmetric relatively the symmetry plane. To 
select the first (symmetric in x — )■ —x transformation) type, we take Cmn = (— l)™'Cmn- As 
a result, we shall obtain the system of equations: 



p- m ^^^^Yl E ^^lAl^ = 0- (20) 

To select the second (antisymmetric in a; — —x transformation) type of plasmon modes 
in Eq. (fT9|) . we take Cmn = — {—l)™'CmL As a result, we shall have 



-ie-i) p-: (tu) i^p^Y, E = o. (21) 

q=0 p=0 

In Fig. m the dependence of normalized plasmon frequencies u/upi = Ij \J\ — e of a 
cluster of two identical oblate nanospheroids on normalized distances //2a between the 
nanoparticles' centers, is shown for the first three plasmon modes. Eigenvalues e were 
obtained as a solution of the equations systems (!20|) and (|2T|) . 
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One can see in Fig. H] that in clusters of two oblate nanospheroids, modes of "T" , "M" , 
and "L" types, that are analogous to T-, M- and L- modes of a cluster made of two prolate 
spheroids (see Fig. [2]), can exist. The T- and M-modes are the solutions of the system (12U1) . 
while the L-modes are the solution of the system fl^T]) . T- and L-modes can be derived 
by the method of hybridization of plasmonic modes of two oblate nanospheroids and their 
plasmonic frequencies are lying in the range < u < Upi / \pl. An infinite number of 
plasmonic frequencies of higher L-and T- modes lie near ujpi/\/2. When width of the gap 
decreases to 0, the ratio oj/ujpi for T-modes tends to various values in the range from 
to 1/ v^, while plasmonic frequencies of L-modes approach zero by analogy with L-modes 
in a cluster of two spherical nanoparticles-^^'^^. Plasmonic frequencies of strongly localized 
M-modes (Fig. Hl^a)) lie in the range ujpi/\/2 < to < cupi, as it happens in a cluster of two 
spherical nanoparticles^^i^. As width of the gap between oblate nanospheroids decreases 
to zero, plasmon frequencies of M-modes tend to bulk plasmon frequency cOpi analogously 
to the case of a two-sphere cluster—*^. For large distances between the spheroids, M- 
modes disappear, and plasmon frequencies of L- and T-modes of a cluster of two oblate 
nanospheroids tend to plasmonic frequencies of a single spheroid (see. Fig. Hl^c)) and can be 
found by means of a self-consistent model with approximation of spheroids by anisotropic 
point dipoles. 

In Fig. O the distribution of a surface charge of plasmon modes of lower orders in a 
cluster of two identical oblate nanospheroids is shown. One can see in this figure that the 
charge distribution is symmetric in T- and M-modes, while in L-mode it is antisymmetric in 
agreement with the definition of these modes. For small distances between nanospheroids, 
charges in M- and L-modes are strongly localized near the gap. On the contrary, when 
the distance between the spheroids increases, the charge distribution tends to symmetric or 
antisymmetric combination of a surface charge in a single oblate nanospheroid (see Fig.[5]^c)). 

Thus, in a cluster of two oblate or prolate spheroidal nanoparticles, fundamental sym- 
metric and antisymmetric plasmon modes of "T" , "M" , and "L" types can be excited, and 
it is these modes that define all optical properties of a two- nanospheroid cluster. 
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III. A CLUSTER OF TWO METAL NANOSPHEROIDS IN THE FIELD OF A 
PLANE ELECTROMAGNETIC WAVE 



In this section, we will consider a two-spheroid cluster in a uniform electric field with the 
potential 



where the time factor e~"^* is omitted. This case corresponds to a plane wave incidence 
and is important for transformation of far fields into near fields, for enhancement of electric 
fields, and for effective excitation of atoms and molecules. Here again, we restrict ourselves 
to two special cases of nanospheroids' shape and position which are most interesting for 
applications. In the first case, two prolate nanospheroids have a general axis of rotation z, 
and axes x and y are parallel (Fig. [T](a)). In the second case, two oblate nanospheroids have 
parallel axes of rotation z and general axis x, i.e., the centers of the nanoparticles are located 
in the plane that is perpendicular to the axes of rotation (Fig. [Hb)). The distance between 
the nanospheroids' centers is more than a half-sum of the lengths between their foci (the 
nanoparticles are not overlapping). 

A. A cluster of two prolate nanospheroids 

Here, we also will use local systems of spheroidal coordinates, the origins of which 
are placed in the nanospheroids' centers (see Fig. [H^a)). The potential inside the j-th 
nanospheroid again can be presented as a series in spheroidal harmonics (Q, while the 
potential outside the nanospheroids now should be presented in the form: 



where v'""*, <f2^^ are contributions from the first and second nanospheroids (see Eq. (ITU]) ), 
and (fo is the potential of the external electric field fl22|l . 

Electric potential of the incident plane wave ( 122|) in local coordinates of the j-th (j = 1, 
2) prolate nanospheroids, looks like: 



(fo = -Eo^x - Eoyy - EozZ, 



(22) 



,out 




(23) 
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^'i^ =f {Eo, COS + Eoy sin (</>,)) (0) (Vj) 

-fE,,P, (0) Pi iVj) + i?o4' (24) 

To find electric potentials and electric fields inside and outside the spheroids, one should 
use the boundary conditions for the potential: 



( , ^out I , „out I 



56 



6=Co 
6=€o 

6=?o 



56 



,„out I ,„out I ,-(2) 



56 



1=«0 



(25) 



where 6 = c/ vc^ 



c/ / is local radial coordinate that defines surfaces of the nanoparti- 



cles, e is permittivity of the nanoparticle's material. Making use of the boundary conditions 
(!25|) and the addition translation theorem (fT2l) . one can obtain the following system of equa- 
tions for the coefficients Cmn , Cmn, Dml. and Dmli that define the outside field (see Eq. (fTOjl 
and Eq. ^) (n = 0, 1, 2, . . . ; m = 0, 1, 2, . . . , n): 



dP^ {^q) pm /£■ n (6) \ ^(1) 

— ^6 — ~ — ^6 — / 



+ (£ — 1) (,^o) S^mnmq^^q— ^^mn ~ 1) (6) 

ct^o — "^0 

q=m 

■ -77. Vn " Uoj ^77 I (-^mn 

"■50 "^02 / 



-J_Cc- 1 A P'Tt /-(^ A (6) c'(l) ^(1)_ (2)/ -I \ pm \ ^^^^^6) 



d^o d^o 



and (n = 1, 2, 3, . . . ; m = 1, 2, 3, . . . , n) 



(26) 
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dPn (^o) ^m u \ pmtc \ i^o) \ n(l) 
^7 UoJ - UoiJ ^7 ^mn 

+ - 1) (eo) ^^^^^E 5(?L.^S= e - 1) (^o) 

+ ie-l) PT (eo) ^^^E CU^S= e ie - 1) (eo) 

5=m 



(27) 



where the coefficients 



= b^l = -5mi6nifEoy, (28) 

define excitation field. In Eq. 0281) . 5ni is Kronecker delta symbol. It should be noted here 
that due to axial symmetry of the considered cluster, the systems of equations fl26|) and fl27|) 
allow one to find the coefficients Cmn and -Dmn for given order m, while degree n runs over 
n = m, m + 1, m + 2, . . . , m + N, where is a large number that defines accuracy of the 
solution. 

Induced dipole moment of a cluster of two prolate nanospheroids, placed in the field 
of a plane electromagnetic wave, can be calculated by analogy with a single prolate 
nanospheroid^i, that is, by finding far field asymptotes of the potential (!23l) . As a re- 
sult, for dipole moment induced in the j-th nanospheroid {j = 1, 2) we have 

4^) = ^CK\ 4^) = ^/^;?, d^^ = ^Ci^, (29) 



and the total dipole moment of the cluster will be the sum of the momenta (12911 . Scattering 
and absorption cross-sections can be easily found if the dipole momenta (129|) are known^: 

,..^,j-yj^^ ,.».=?^f-yE;, (30) 

\vj |E„|' 3 \vj \eX 

where d = d^^^ + d*-^) denotes dipole momentum of the whole system, and asterix denotes 
an operation of the complex conjugation. 
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In Fig. [6l absorption (a) and scattering (b) cross-sections of a cluster of two identical 
prolate nanospheroids made from silver, are shown as a function of the wavelength. For 
longitudinal (z) polarization, both of the cross-sections have two peaks that correspond to 
longitudinal plasmonic oscillations with L = 1, 2 (see Fig. [2]^b)). It is very important that 
both of the peaks are split substantially relatively the case of a single spheroid (the "z" 
dashed curve) due to strong interaction between the nanospheroids. 

On the contrary, for transversal (x or y) polarization one can see only one peak due to 
excitation of the symmetrical T = 1 mode, and this peak is shifted just slightly relatively 
the single spheroid resonance (the "y" dashed curve). It means that transversal (x or y 
polarization) excitation of a two-spheroid cluster induces only weak interaction between 
the nanospheroids (see the dispersion curves for T-modes on Fig. El^a)). Due to this weak 
interaction, absorption and scattering cross-sections are approximately equal to doubled and 
quadrupled cross-section of a single spheroid respectively. 

It should be noticed that in Fig. |6] the maxima of absorption, corresponding to plasmon 
oscillations of M-type that should he in the interval ijjpi/\/2 < u < Upi, that corresponds 
to 326 < A < 337 nm for silver^, are not visible. It is related to the fact that M-modes 
interact with a homogeneous electric field weakly and can be effectively excited only by a 
source of radiation that is nonuniform in comparison with a size of the gap between the 
nanoparticles^ (see Fig. |3]). 

To control correctness and accuracy of our analytical calculations, we have also carried 
out finite element simulation of this system with Comsol Multiphysics@ software. The 
results of this simulation are shown by the circles in Fig. O One can see that there is a 
fine agreement between the analytical and pure numerical calculations. This fact confirms 
correctness and accuracy of the both of approaches. 

B. A cluster of two oblate nanospheroids 

The case of two oblate nanospheroids is in many aspects similar to the case of two prolate 
nanospheroids, considered above. So let us again choose local systems of coordinates that are 
connected with each of the nanospheroids and have origins in their centers (see Fig. lU^b)). 
The potential inside the j-th nanospheroid again can be presented as a series in spheroidal 
harmonics ( 1T6|1 . while the potential outside the oblate spheroids can be presented in the 
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form: 



,-OUt ,^out I ,„OUt I ^Q1^ 

f =^1+^2+ fo, (31) 
where y9°"*,v^2"* contributions from the first and second nanospheroids (see Eq. f lT7|) ) 
and (fiQ is the potential of the external electric field (12^ . In the local coordinates of the j-th 
(j = 1, 2) oblate nanospheroids, it looks like 



y.?) =-tf {Eo, cos (0,) + Eoy sin (0,)) P^^ {zQ (r^,) 

+z/Eo.Pi (^e,) ^1 iVj) + (-1)' ^o4- (32) 



Making use of the boundary conditions (1251) with = c/ \/ — = c/ f and the addition- 
translation theorem (!T8|l . we shall obtain the following systems of equations (n = 0, 1, 2, 



m = 0, 1, 2, . . . , n) 



+ (. - 1) (z^o) ^^^f^ E ^i'4.^i?= - 1) Pn (^eo) 

+ (. - 1) P- (z^o) y E M^LciS= ^ - 1) Pn i^^o) 



and (n = 1, 2, 3, . . . ; m = 1, 2, 3, . . . , n) 



+ (^ - 1) i^^o) ^^^E E ^a.4?= e is - 1) pr (.^o) 
+ (^ - 1) Pn i^^o) ^^^E E ^i^.^S^= e - 1) (^eo) 



5 



(34) 
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where 



€i = €i = ^SmlSnlfEoy. (35) 

Apparently, the equations (!33|) and (134|) have more comphcated structure than Eqs. (!26|) 
and fl27|) because now, due to lack of axial symmetry, one can not split the system of 
equations into systems with fixed order m of the Legendre function. 

For calculation of absorption and scattering cross-sections of a cluster in the field of a 
plane electromagnetic wave, one can again use Eq. fl30|) . where dipole momenta of each 
spheroid can be expressed through the solutions of Eqs. (l33l) and flMl) (j = 1, 2) 

4^) = -^cK\ 4^) = -^/^;?, d^^ = -^Ci^. (36) 

In Fig. [71 the absorption (a) and scattering (b) cross-sections of a cluster of two iden- 
tical oblate nanospheroids made from silver are shown as a function of wavelength. For 
longitudinal (x) polarization, both the cross-sections have two peaks, that correspond to 
antisymmetric plasmonic oscillations with L = 1, 2 (see Fig. HJ^b)). It is very important that 
now only one peak (L = 1) is shifted substantially relatively the case of a single spheroid 
(the "x, y" dashed curves) due to strong interaction between the nanospheroids. The L = 
2 mode suffers only a small shift in agreement with Fig. Hl^b). 

For transversal (y) polarization, one can see only one peak due to excitation of the sym- 
metrical T =1 mode, and this peak is only slightly shifted relatively the single spheroid 
resonance (the dashed curve). It means that transversal (y polarization) excitation of a 
two-spheroid cluster results only in weak interaction between the nanospheroids (see the 
dispersion curves for the T- modes on Fig. 111(a)). Due to this weak interaction, absorption 
and scattering cross-sections for this polarization are approximately equal to doubled and 
quadrupled cross-sections of a single spheroid respectively. It is also interesting that plas- 
monic frequency of the L = 2 mode is very close to plasmonic frequency of the T = 1 mode. 
This fact can be easily understood from analysis of Fig. |H Indeed, when width of the gap 
tends to zero, plasmonic frequency of L = 2 modes also decreases to zero, while plasmonic 
frequency of the T = 1 mode increases slightly. So, at some point these modes will intersect 
and have the same frequencies, and we observe this situation in Fig. [3 
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It should be noticed that in Fig. [7] the maxima of absorption corresponding to plasmon 
oscillations of M-type again are not visible. It is related to the fact that M-modes interact 
with a homogeneous electric field weakly and can be effectively excited only by a source of 
radiation that is nonuniform in comparison with a size of the gap between nanoparticles^ 
(see Fig.E]). 

To control correctness and accuracy of our analytical calculations for a cluster of two 
oblate spheroids, we have also carried out finite element simulation of this system with 
Comsol Multiphysics® software. The results of this simulation are shown by the circles in 
Fig. [71 One can see that there is a fine agreement between the analytical and pure numerical 
calculations. This fact confirms correctness and accuracy of the both of approaches again. 

C. Enhancement of local fields 

The most important characteristic of nanoparticles clusters is the factor of incident field 
enhancement in the gap between the nanoparticles. It is the characteristic that allows 
determining the excitation rate of molecules near the nanoparticles or intensity of surface- 
enhanced Raman scattering-. Moreover, achievement of high values of this factor is the 
main goal of optical nano- antennas development. 

Distribution of squared electric field for the L=l resonance in a cluster of two prolate 
spheroids is shown in Fig. [HI that shows that, indeed, the maximal field enhancement takes 
place in the gap between the nanoparticles on their surfaces. The field maxima are also 
present on the outer side of the cluster; however, field amplitude is essentially less there. 
According to general theorems for harmonic functions, the field maximum can be reached 
only on the region boundaries. In our case, the field maxima are reached in those points of 
the spheroids' surface where the distance between the spheroids is minimal. 

Using Eqs. ( ITOl) and ( IT71) . one can find explicit expressions for the field enhancement 
factor G. For clusters of two identitcal prolate spheroids in the considered configuration 
(Fig. [H^a)), one can obtain the following expression for the field maximum in the case of an 
incident field polarized along the z-axis: 
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1)" Cq^ . In the most interesting 



where ii = c/ f and ^2 = — c) / /, while Co„ 
case of a small gap and strongly prolate spheroids, ^1, ^2 ~ 1, and one may use the asymptotic 



form 



2(5-1)' 



As a result, the field enhancement factor takes the form: 
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(38) 



In the case of clusters of two identical oblate spheroids (Fig. [U^b)) and incident field 
polarized along the x-axis, we obtain 
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(39) 



I J- ) '^mn — '^mn- 



where 6 = V«V/' " 1 6 = ^(/-a)^/^-!; C^m,! 

In Fig.[9l the dependence of squared electric field enhancement (!37|) and (!39l) for clusters of 
two identical silver nanospheroids on the wavelength, is shown. Comparing peaks positions 
to the dispersion curves in Fig. [2] and Fig. HI one can come to a conclusion that only "L" 
types of plasmon modes are excited in the clusters for the considered configurations of 
nanospheroids and incident electromagnetic wave polarizations (along the line joining the 
nanoparticles' centers). In particular, excitation of the L = 1 and L = 2 modes is noticeable. 
At that, the position of squared field enhancement peaks agrees with the maxima of the 
absorption and scattering cross-sections shown in Fig. [H] and [7] by the solid lines z and x 
correspondingly. It should be mentioned that the value of squared field enhancement near 
a cluster of two nanospheroids can reach up to 10®. In the case of single nanoparticles, 
this value is almost two orders less than that of clusters (cf. the solid and dashed curves 
in Fig. [9]). This fact determines a greater attractiveness of metal nanoparticles clusters in 
comparison to single nanoparticles for investigation of Surface Enhance Raman Scattering 
(SERS) and Surface- Enhanced Fluorescence (SEE). Let us note that the obtained great 
field enhancement factors can be slightly less in practice since for small particles and for 
small gaps between them, non-local and other effects not considered in this research become 
essential. 
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IV. A CLUSTER OF TWO NANOSPHEROIDS IN THE FIELD OF A DIPOLE 
SOURCE 



In the previous section, we have considered the case of a nano-antenna placed in the 
field of a plane wave. However, highly nonuniform optical fields occur very often in nano- 
environment. For example, such fields arise when a plasmonic nano-antenna is excited by 
an atom, or a molecule, or another nanolocalized source of light. So, in this section, we 
will consider an important case of a two-nanospheroid cluster in the field of electric dipole 
sources. Excitation of the cluster by magnetic dipole and electric quadrupole sources can 
be analyzed analogously. 



A. A cluster of two prolate nanospheroids 

A case of two prolate nanospheroids in the field of a dipole source of radiation can be 
considered by perfect analogy with a case of the same cluster in a uniform field. One should 
again look for solution in the form (Q, (fTOj) and then apply the boundary conditions (!25|l . 
The only difference is that now the external potential is the potential ipo of the dipole that 
has the following form in the j-th local system of coordinates of a prolate spheroid^ {j = 1, 
2): 




n=0 m=0 



1 

« (0) Pn {m) 



(doV0 a^^n cos (m0j) + (doV0 sin (m0 
(doV^) 7mn COS + (doV0 (Tmn siu (m0. 
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In Eq. (HOj) . do denotes the dipole momentum of a source placed at r', is a gradient over 
r' in local coordinates, and 
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are expansion coefficients of the unit charge potential in local coordinates of a prolate 
spheroid. 

As a result of applying of the boundary conditions, one can obtain a system of equations 
for the unknown coefficients Cmn, Dmn in Eq. ( fTOl) . The new system can be easily derived 
from Eqs. fl2^ and fl7r|) if one makes the following replacement for the coefficients amn and 



h^l = - (doV;) /3« , = - (doV'2) /^i^) , (42) 

where amn and /3mn are defined by Eq. ( |4T|) . 

After Eqs. (126|) and (!27|) have been solved with taking (142!) into account, one can find the 
total induced dipole moment of the prolate nanospheroids 



dx — {c[i + C[i 



p 

d, = — 
3 



{Cl\^+Ci?). (43) 



Knowing the dipole momenta fl43l) . it is easy to find radiative decay rate of a dipole placed 
near the prolate nanospheroid^: 

prad , .3 

where P^"''^ is radiation power at frequency u and hu is emitted photon energy. 

Radiative decay rate is a very important characteristics in such applications as surface 
enhanced Raman scattering, surface-enhanced fluorescence, nanolasers, and so on. To char- 
acterize the radiative decay rate, it is naturally to normalize it to radiative decay rate of a 
dipole in free space, 70 = ^ = ^\dof. 

In Fig. [ini normalized radiative decay rate of a dipole source placed at the middle point 
of the gap is shown. As it is well seen in Fig. fTOT a). if the distance between the prolate 
nanospheroids is small (Fig. [TOT a). the curves a and 6), the dipole source with a moment 
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oriented perpendicular to the cluster's axis of rotation, can excite both symmetrical T- and 
M-modes. This fact differs from the case of excitation of the same cluster with a plane 
wave, when M-modes with peaks located in the region of A < 337 nm (see Fig. [2]^a)) are 
not excited. When the distance between the nanospheroids increases (see Fig. [TUT a)). the 
peak corresponding to M-modes shifts to A ~ 337 nm (w ~ ujpi/\/2) and then disappear. 
After that point, only the peaks corresponding to plasmonic T-modes can be observed. Of 
course, this picture is in agreement with the behavior of the plasmonic M-modes shown in 
Fig. [2t^a). We shall also notice that for large enough distances between the nanospheroids 
(see Fig. [TUlfa). the curve 7) the self-consistent model^^i^, in which nanoparticles are replaced 
by point dipoles with corresponding polarizabilities^*^, can be effectively used for calculation 
of radiative decay rate (the dashed curve). 

When dipole moment of a source is oriented along the axis of symmetry (Fig. [TUT b)). 
only antisymmetric L-modes can be excited due to symmetry reasons. From Fig. [TOT b). one 
can also see that for small enough distances between the spheroids there are two plasmonic 
modes (L = 1, 2) that interact with the dipole source. When the distance between the 
spheroids diminishes, further peaks of radiation power shift towards long wavelengths. At 
large distances between the nanospheroids, there is only one maximum corresponding to the 
L = 1 plasmonic mode (see Fig. [TOT b). the curve 7). In this case, the radiative decay rate 
of a dipole placed near a cluster of two prolate nanospheroids can be also calculated with 
making use of the self-consistent analytical model in which the spheroids are approximated 
by point dipoles (see the dashed curve in Fig. [TOT b)). 



B. A cluster of two oblate nanospheroids 

A case of two oblate nanospheroids in the field of a dipole source of radiation can be 
considered by perfect analogy with a case of the same cluster in a uniform field. One should 
again look for solution in the forms ( fT6l) and ( ITTl) and then apply the boundary conditions 
( I25II . The only difference is that now the external potential is the potential ip^ of the dipole 
that has the following form in the j-th local system of coordinates of an oblate spheroid^ 
{j = 1, 2): 
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In Eq. (HSj) . do denotes the dipole momentum of a source placed at r', V^- is gradient over 
r' in local coordinates, and 
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are expansion coefficients of the unit charge potential in local coordinates of an oblate 
spheroid. 

As a result of applying of the boundary conditions, one can obtain the system of equations 
for the unknown coefficients Cml, Dmn in Eq. f|T7j) . The system of equations for these coeffi- 
cients can be easily derived from Eqs. (!33|) and (!34l) if one makes the following replacement 
for the coefficients amn and hn 



Jmn 



a. 



mn = -i (doV'i) a^i, a^^^ = -i (doV'ai) a^^i 



feW = -I (doV;) f3^2, b^l = W,) (3i^l (47) 

where amn and /3mn are defined in Eq. (B6l) . 

The total induced dipole moment of a cluster of two oblate nanospheroids can be found 
by using the following expressions (cf. Eq. f HSj) ): 
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In Fig. [m normalized radiative decay rate of a dipole placed at the middle point between 
the oblate nanospheroids is shown. As one can see in Fig. fTTT a). if dipole moment of a source 
is oriented along the line connecting the oblate nanospheroids' centers, only plasmonic L- 
modes are excited as it took place for the case of a cluster of two prolate nanospheroids (cf. 
Fig. [TUT b)). When the distance between the nanospheroids is large enough (the 7 curve), 
the properties of radiative decay of a dipole source located near the cluster can be more 
and more precisely approximated by means of the self-consistent model where spheroids are 
modeled by point dipoles (see Fig. [TTT a). the dashed curve). 

From Fig. [TTT b) and (c) that correspond to the case of dipole moment of a source oriented 
along the y and z axes, one can conclude that only symmetric T-modes are presented here. 
Nevertheless, at small enough distances between the nanospheroids one can expect that 
M-modes are also excited but can not be seen due to large losses (large imaginary part of 
dielectric permittivity) in silver spheroids. If the imaginary part of permittivity of spheroidal 
nanoparticles' material is small enough, for example, in a case of silicized carbon (SiC)^, 
all fine features of the spectra and M-modes in particular will be well visible as it was 
demonstrated in2S>^ by the example of a cluster of two spheres. 

V. CONCLUSION 

Thus, in the present work optical properties of clusters made of two metal nanospheroids 
are considered theoretically, and analytical results are obtained. This investigation has 
become possible due to proving of the addition translation theorem for spheroidal functions 
in quasistatic regime. Plasmonic eigenoscillations were analyzed in details, and it was found 
that in a cluster of two prolate or oblate nanospheroids there can be three types of plasmon 
modes. Two of them (low frequency, < a; < ujpi/\/2, L- and T-modes) can be effectively 
excited by a plane electromagnetic wave, while the third type (high frequency, ijjpi/\/2 < tu < 
Upi, M-modes) can be excited only by a strongly nonuniform field of a nano localized source 
of light (a molecule, a quantum dot) located in the gap between two adjacent nanoparticles. 

We have also investigated excitation of a nano-antenna made from two silver nanospheroids 
by fields of a plane wave and an electric dipole. The results of these investigations allow 
us to find absorption and scattering cross-sections of the nano-antenna as a function of 
wavelength for various polarization of an incident plane electromagnetic wave and to at- 
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tribute all observable peaks to excitation of corresponding plasmonic modes. We have also 
analyzed radiative decay rate (or local density of state) of a dipole placed in the gap between 
the nanospheroids and have attributed all observable peaks to excitation of corresponding 
plasmonic modes. 

The obtained analytical results can be used in many applications based on plasmonic 
nano-antennas and on enhancement of local field (SERS, SEP, nanolasers, nano-optical 
circuits, and so on). Besides, our results are very important for controlling accuracy of 
different computational software that has no a priori test of accuracy. 
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Appendix: Translational addition theorem for spheroidal wave functions in a qua- 
sistatic Umit 

A solution of the Helmholtz equation in spheroidal coordinates can be expanded as a 
series in spheroidal wave functions^. In spite of the fact that these functions are thor- 
oughly studied^""— 1^ , the problem of electromagnetic waves scattering on a spheroid still 
remains one of the most complex ones, first of all, due to mathematical aspects related to 
usage of spheroidal wave functions. On the other hand, it is well-known that in a case of 
spheroids which size is substantially less than the wavelength (quasistatic limit), a solution 
of the Laplace equation in spheroidal coordinates can be presented as a series in associated 
Legendre functions^-. The Legendre functions are widespread and included as standard spe- 
cial functions in majority of mathematical software programs, such as Maple, Mathematica, 
MatLAB, etc. That is why it is important to obtain the translational addition theorem for 
spheroidal wave functions in the case when sizes of spheroids are far less than the wavelength 
in which these functions would be "replaced" by corresponding Legendre functions. 

Iri^, a rotational-translational addition theorem was obtained for spheroidal wave func- 
tions in a case of an arbitrary size of spheroids in comparison to the wavelength. This 
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theorem was derived in the following manner: first, spheroidal wave functions were ex- 
panded over spherical wave functions, then the spherical wave functions were expanded over 
spheroidal ones. Substituting one expansion into the other and using the addition theorem 
for spherical wave functions^i^, the required theorem for spheroidal functions was derived 
as a result. Here, we will do the same in a particular case when spheroids' sizes are far less 
than the wavelength. 

In of a prolate spheroid, we have the following expansions^ 
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- ^ '^''-''d^-\m\W)h';^Hkr)P^{cos9)e^-t (A.l) 

q=\m\,\m\+l 

where 1 < ^ < oo, —1 < 1] < 1 and < < 27r are coordinates of the prolate spheroid^; 
< r < oo, < < TT and < < 27r are spherical coordinates; Rmn{kf,^) and 
Sinn {kf, rj) are radial and angular spheroidal functions of the third and first kinds according 
to^; (kr), and (cos^) are the spherical Hankel function of the first kind^^ and the 
associated Legendre function-S; c?^^| (kf) are the expansion coefficients^; k is wavenumber, 
/ is a half of the spheroid's focal distance; the stroke of the sum character implies summation 
over only even or only odd values of q. Another important relation from^ has the form: 



j„ (kr) P- (cos 6) e-^ = 

2n + l[n — m)\ 



ij=|m|,|m|+l 

where jn (kr) is the spherical Bessel function^; Rmg (fc/, is a radial spheroidal function of 
the first kincl^, and Nmq {kf) is normalization factor of angular spheroidal functions of the 
first kind^. In the case of a spheroid which size is far less than the wavelength, that is, in 
the case /c — )■ 0, one can use the following asymptotic expressions^: 
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and Nran ~ 2^ (n-m)! " Radial and angular spheroidal functions Rmn{kf,C)-, Rmn{kf,C.) 
and Smn {kf, rj) can be also presented in the asymptotic fornt^i^ 



Rl^n (kf, ~ (2^ -l)!!(2r! + l)!! ' 

( v. 'rt « -C (>)) . (A.4) 

where (77) is the associated Legendre function^ defined on the segment —1<1]<1] 
-PjT (O5 Qri (0 are associated Legendre functions of the first and second kinds^, correspond- 
ingly, defined in a complex plane with branch cut from —00 to +1. Finally, asymptotic forms 
for spherical Bessel and Henkel functions have the form^: 

Jn M ^ -^^^ (fcr)" , /i« (kr) ^ -ziM: (A;r)-"~^ . (A.5) 
JnK J (2n+ 1)! ^ ^ n K J 2"^! ^ ^ ^ ^ 



Making use of Eqs. (1A.3I) - (1A.5I) . one can obtain the following relations from Eqs. (lA.lll 
and dQ): 
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To obtain the required addition theorem for wave functions of prolate spheroids which 
size is far less than the wavelength from Eq. ( ]A.6p . it is necessary to apply the corresponding 
theorem for spherical wave functions^: 
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In the expression f lA.7p . the indices j, s = 1, 2 (j 7^ s) denote the first "1" and the second "2" 
local systems of spherical coordinates, that are based on local Cartesian coordinates with 
parallel axes and attached to the spheres' centers; (r^s, Ojg, (pjg) are spherical coordinates of 
the origin of the s-th coordinate system in the j-th local system of coordinates; / = vjg = r^j 
is the distance between the origins of local coordinates systems (see Fig. [T^ . Substituting 
Eq. ( 1A.7P into Eq. (lA.6p . we obtain after series of transformations the required translational 
addition theorem for wave functions of prolate spheroids in the quasistatic limit: 



00 q 



9=0 p=-q 

{-If {2q + I) {q - p)\ {n + m)\ 



00 00 



{q + p) \ {n — m)\ 
{q + n + 2r + 2k — p + m)\ 



r=o fc=o + 2^ + 1)" (2^ + 2A; + 1)!! (2r)!! {2k)\\ 

2r+q / /. \ 2k+n+l 

I .Ik \ 

X 



j) [f) n+n";2.+2. (cos e,,) e^(--)^- . (A.8) 

In Eq. (lA.Sp . the indices j, s = 1, 2 (j 7^ s) denote the first and second local prolate spheroid 
coordinate systems with the parallel Cartesian axes introduced above; the spherical angles 



9js and ipjs are the same as in Eq. ( ]A.7I) : I is the distance between the origins of local 
coordinate systems (see Fig. [12] for details). Let us add that an obligatory condition for 
usage of the theorem flA.SP is parallelism of Cartesian axes in the local coordinate systems. 



In the case when the second local coordinate system is turned around the first one, the 
theorem flA.7l) is to be changeci^, resulting in corresponding modification of the theorem 



(1A.8I) . Let us also note that mathematical applicability condition for ( lA.SI) is convergence 
of the series in this expression. In particular, one may state that Eq. ( lA.Sp is applicable in 
the case of non-overlapping spheroids^. 

The translational addition theorem for functions of oblate spheroids which size is much 
less than the wavelength, can be obtained from Eq. ( lA.SI) if the formal substitution is usecl^: 
^ — )■ and / — )■ —if. In expressions derived in this manner, ^ is a coordinate in the system 
of the oblate spheroid (0 < < 00); / is a half of the focal distance of the oblate spheroid. 
Let us note that (1A.8I) is a very general theorem, and many other theorems relating prolate, 
oblate, and spherical functions can be derived from it. 
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The obtained theorem (lA.SP is not very convenient for numerical calculations since it 
contains summation both over positive and negative p. To simplify it, let us use properties 
of the associated Legendre functions^i^, basing on which we shall find that 



■1) 



p—m 



{q + p) \ {n — m)\ 



1 2 



{q — p)\{n + m)! 

Combining Eqs. flA.Sp and flA.9p . one can derive the following relation (m, p = 0, 1, 2, . . . ): 



{cos imcbj) 
sm(m0j) g=op=o 

) COS {p(f)s - (p-m) (pjs) 



X 



Lpqmn {fj, fs, I, Ojs) siu {p(j)s - {p - m) (j)js) 

+Tpqmn {fj, fs, I, Ojs) COS {-pcps + {p + m) (pjs) 
+Tpqmn {fj, fs, I, Ojs) siu (-^0^ + {p + m) (f)js) 



(A. 10) 



where 



, . ^ {-!)' {2 - 6op) {2q + 1) {q - p)\ {n + m)\ 

J-'pqmn \Jji Jsi "jsj 



2{q + p) \ {n — m)\ 
[q + n + 2r + 2k — p + m)! 
^=0 fc=o + 2r + 1)!! {2n + 2k + 1)!! (2r)!! (2A;)!! 

2r+q / _f \ 2k+n+l 



oo oo 

X 



^ , (-ir+-(2-M(2g + l)(g-p)!(r^ + m)! 

-I- pqmn \ Jji J si ''js/ 



oo oo 



2 (g + p)! — m)! 
(g + n + 2r + 2A; — p — m) ! 



,.=0 fe-o + 2^ + 1)" (2"- + 2A; + 1)!! (2r)!! (2A;)!! 



yj (fj P^:Z2r^2,i^osO,s), (A.ll) 

and ^op is Kronecker delta symbol. 

In the case of coaxial prolate spheroids with a common axis z, one should take = (ps 
and Ojs = or Ojs = tt in Eq. f lA.10|) . This results in the following relation (m = 0, 1, 2, 
...): 
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Qn (0) PT {m) = E (/^' ^' ^^•^) C (^^) n'" (^«)' 

q=m 

-iy{2q + l) (g-m)!(n + m)! 
(g + m)! (n — m)! 

(g + n + 2r + 2fc)! ( U\ f fA 

(2g + 2r + l)!!(2n + 2A; + l)!!(2r)!!(2A;)!! V / / \l) 



iJmqmn \Jji Jsi ^js) /' i M /' M ^q+n l^COS t/j^J 



(A.12) 



In the case when local coordinate systems of the spheroids have a common axis x, one 
should take 9js = 7r/2 and (j)js = or (f)js = tt in Eq. (lA.lOp . As a result, we have in this 
case 



Qn iQ Pn im) 



COS [m(pjj 



sm 171(1 



I -^pqmn {fjy fsi 

I, (pjs) sin' 



q=0 p=0 



where 



To simplify Eq. ( ]A.14p . it is convenient to use the relation^: 



(A.13) 



-^pqmn ifj: fsi ^i 'Pjs) 

ifji fsi li 71-/2)) COS ((p - m) (pj,) , 

-l^pqmn ifji fsi ^i 'Pjs) 

= {Lpgmn ifji fsi li 71-/2) - Tpqmn (fji fsi U Vr/2)) COS ((p - m) (f)js) ■ (A.14) 



^ ^ 2™ /, , 7r\ r((n + m + l)/2) , 
Pr(0) = ^coB((„ + ™)-)^|L- (A.15) 

where F (x) is Euler's gamma function^. 
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FIG. 2. (Color online) Normalized frequencies of the first three plasmonic oscillations in a clus- 
ter of two identical prolate spheroidal nanoparticles as a function of the distance between the 
nanospheroids' centers. The axis-symmetrical case (m = 0) is considered, and the aspect ratio of a 
single spheroid is a/c = 0.6 (^o = l/y^l — (a/c)^). (a) Symmetrical modes (eigenvalues of Eq. (fT4|) ). 
(b) antisymmetric modes (eigenvalues of Eq. (jlSp ). Dashed lines show plasmon frequencies of a 
single prolate nanospheroid ^ = {^^g^Q^T (^o) / {^^^^Qn (^o) - (Co) ^^|^) }'^'- The 
figure (c) shows plasmon frequencies of the single prolate nanospheroid as a function of the inverse 
aspect ratio c/a. The vertical line corresponds to a/c = 0.6 and allows to select asymptotic values 
for the figures (a) and (b). 
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FIG. 3. (Color online) Distribution of a surface charge (a.u.) of the lowest plasmon mode in 
a cluster of two identical prolate spheroidal nanoparticles according to solution of the equations 
(jl4p and (jl5p . The axial-symmetrical case {m = 0) is considered, and the aspect ratio of a 
single spheroid is a/c = 0.6. The distance between the nanospheroids' centers is //2c = 1.03 (a), 
//2c = 1.2 (b). Distribution of a surface charge in a single spheroid is shown in panel (c). The red 
color corresponds to the positive charge, and blue - to the negative one. 
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FIG. 4. (Color online) Normalized plasmon frequencies of the first three modes in a clus- 
ter of two identical oblate nanospheroids as a function of the distance between their centers, 
(a) Symmetrical modes (the eigenvalues of Eq. (|20p). (b) antisymmetric modes (the eigenval- 
ues of Eq. (|21|) ). By dashed lines, plasmon frequencies of a single oblate nanospheroid — = 
I d^Mg™ (i^o) / (^^Em^^Q^ (i^o) _ (i^o) ^^Q^) are shown. The spect ratio of a 
single oblate nanospheroids is c/a = 0.6 (^0 = ~ The figure (c) shows plasmon fre- 

quencies of a single oblate nanospheroid as a function of the inverse aspect ratio a/c. The vertical 
line corresponds to c/a = 0.6 and allows to select asymptotic values for the figures (a) and (b). 
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FIG. 5. (Color online) Distribution of a surface charge (a.u.) of the lowest plasmon modes in 
a cluster of two identical oblate spheroidal nanoparticles according to solution of the equations 
(j20p and (|2ip . The aspect ratio is c/a = 0.6. The distance between the nanospheroids' centers 
I /2a = 1.05 (a), I /2a = 1.4 (b). Distribution of a surface charge in a single spheroid is shown in 
panel (c). The red color corresponds to the positive charge, and blue - to the negative one. 
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FIG. 6. (Color online) Absorption (a) and scattering (b) cross-sections of a cluster of two identical 
prolate nanospheroids made from silver as a function of the wavelength. The large semi-axes of 
the nanospheroids are c = 15 nm, the aspect ratios are a/c = 0.6, and the distance between the 
nanospheroids' centers is l/2c = 1.05. The labels x, y, z correspond to polarization of an incident 
wave along the x, y, z axis. The circles correspond to finite element simulations with Comsol 
Multiphysics(§) software. Doubled absorption cross-section and quadrupled scattering cross-section 
of a single nanospheroid are shown by dashed curves. 
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FIG. 7. (Color online) Absorption (a) and scattering (b) cross-sections of a cluster of two identical 
oblate nanospheroids made from silver as a function of the wavelength. The large semi-axes of 
the nanospheroids are a = 15 nm, the aspect ratios are c/a = 0.6, the distance between the 
nanospheroids' centers is I /2a = 1.05. The labels x, y, z correspond to polarization of an incident 
wave along the x, y, z axis. The circles correspond to finite element simulations with Comsol 
Multiphysics@ software. Doubled absorption cross-section and quadrupled scattering cross-section 
of a single nanospheroid are shown by dashed curves. 
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FIG. 8. (Color online) Spatial distribution of |E|^/|Eor 

in a cluster of two prolate nanospheroids 
for A = 425 nm at longitudinal polarization of the excitation field (L = 1 plasmonic resonance, 
common logarithmic scale). 
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FIG. 9. (Color online) Enhancement of |E|^ in the gap between two identical silver nanospherodis 
as a function of the wavelength. P and O labels correspond to prolate and oblate spheroids 
correspondigly. The large semi-axes of the nanospheroids are 15 nm, the aspect ratios are 0.6, 
I /2c = 1.05 for a cluster of prolate spheroids and l/2a = 1.05 for a cluster of oblate spheroids. 
Enhancements for single nanospheroids are shown by the dashed curve. 
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FIG. 10. (Color online) Normalized radiative decay rate of a dipole placed at the middle point 
between two identical prolate nanospheroids made from silver as a function of the wavelength. The 
dipole source moment is oriented along the x or y axes (a) and along the z axis (b). The large 
semi-axes of the nanospheroids are c = 15 nm, the aspect ratios are a/c = 0.6. The curves a, /3, 7, 
6 correspond to //2c = 1.05, 1.1, 1.3 and 1.03, respectively. The asymptotic expression obtained 
by approximation of the spheroids by point dipoles (//2c = 1.3) is shown by the dashed curve. 
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FIG. 11. (Color online) Normalized radiative decay rate of a dipole placed at the middle point 
between two identical oblate nanospheroids made from silver as a function of the wavelength. The 
dipole source moment is oriented along the x (a), y (b), and z (c) axis. The large semi-axes of 
the nanospheroid are a = 15 nm, the aspect ratio is c/a = 0.6. The curves a, (3, 7 correspond to 
I /2a = 1.05, 1.1, and 1.3, respectively. The asymptotic expression obtained by approximation of 
the spheroids by point dipoles {I /2a = 1.3) is shown by the dashed curve. 
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FIG. 12. Geometry of derivation of the translational addition theorem. 
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